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Abstract 

A rigorous Path Integral construction for a wide class of Weyl 
evolution operators is based on a pseudo-differential O-calculus 
on flat phase spaces of finite and infinite dimensions. 

1 Introduction 

Among the various uses of the Path Integrals the earliest and foremost 
is the application to quantum evolution. The Path Integral formalism 
is appreciated because of its compact operational notation with all the 
gamut of the finite dimensional Integral Calculus, such as integration 
by parts, repeated integration, canonical substitutions, analytic contin- 
uation (Wick rotation), stationary phase approximations etc. 

However, in R.Feynman words, the Path Integral is "an intuitive 
leap at mathematical formalism" . A natural justification would be by 
suitable integral approximations, the route chosen originally by Feyn- 
man himself in the 40's via time-slicing and discretization. Unfortu- 
nately the discretization ambiguities along with the convergence prob- 
lems have plagued the deed from the start (the notable exception |15| 
is for rather special hamiltonians) . 

Nevertheless we propose a rigorous time-slicing construction of the 
(flat) phase space Path Integral for propagators both in Quantum Me- 
chanics and Quantum Field Theory for a fairly general class of quasi- 
dissipative quantum observables (e.g. the Schrodinger hamiltonians 



with smooth scalar potentials of any power growth). Moreover we al- 
low time-dependent hamiltonians and a great variety of discretizations, 
in particular, the standard, Weyl, and normal ones. 

Abstract Cauchy Problem. Consider the Initial Value Problem 

^ + A(t)i(>(t) = 0, m=^o, 0<t<T, #)eH, 

wherein H is a Hilbert space and A(t) is a family of (usually) unbounded 
operators on H. 

The Cauchy Problem is called proper relative a dense subspace S in 
7i if there is the unique solution ip for every ip e S and the Evolution 
Operator 

U{t",t')^{t') = ^{t"), ^(0 e ^, *" > 

is bounded on 7i. 

The Evolution Operator may be sought in the time-slicing form of 
the strong operator limit Product Integral: 



U(t",t') = Yl exp[-A(t)dt] 

t">t>t> 

:= lim U r (t",t') := lim TT exp[-4(t 7 -)Ai 7 -]. 

*j + l>*j 

Here V is a finite partition < £' = £ < • • • < tj < tj+i < • • • < t" — 
t v < T of the interval t' <t<t", Atj = t j+1 -tj, \V\ = maxj |At,-| 



Prom Product to Path Integrals on the Phase Space. Consider 
a quantum evolution equation on £ 2 (IR d ) 

dib % % d 

-jfi(t, q) + -fit, q, i~)^{t, q) = 0, ^(0, q) = ^o, 

with a pseudodifferential operator f(t,q, |^-) in the standard form of 
the qp— quantization of a complex- valued function fit, q, p) on the phase 
space M? d , the standard symbol of f(t, q, |J^)- 

The Tobocman's version of the Dirac-Feynman Ansatz: For small 
At the standard symbol < p\U(t+At, t)\q > of the propagator U(t + At, t) 
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is approximately equal to exp[—^f(t,q,p)At\. Then according to the 

product rule for the standard symbols, the standard symbol of U v is 
approximately equal to the distributional multiple integral 

d-l . d-l 

Y[dX h (tj) exp -^2[p(t j+1 )Aq(tj) - /(^(i^p^+^At,] 

where d\n{tj) = (27rti)~ d dq(tj)dp(tj) is the Lebesgue-Liouville measure 
on the phase space ~R 2d . 

As the mesh \V\ — > 0, the multiple integrals are presumed to con- 
verge to a Hamiltonian Path Integral for the standard symbol of the 
evolution operator U(t", t') 

\{ d\ h (t)exp- / \p(t)q(t)-f(t,q(t),p(t)]dt 

t">t>t> *' 

where d\n{t) = (2irh)~ d dq(t)dp(t) is "the Feynman-Liouville measure" 
on the space of paths from q(t') = q to p(t") = p in M. 2d and 

t" 

[p(t)q(t)- f(t,q(t),p(t)]dt 

is the hamiltonian symplectic action functional on that space. 

The standing physical presumption: all calculus rules are valid in 
the limit. However there are two fundamental mathematical problems: 
The validity of the DFT-ansatz and the existence of the limit. So far 



both problems have been settled only for special / (c.f. [ITR |14], 



Euler Detour. Euler polygonal approximation for the solution of the 
quantum Cauchy Problem 

^ + f{t)m = o 

is the finite difference approximation 

+ f{tj)tl){tj) = 0, V(*o) = ^o, 
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or 

1>(t j+1 ) = (l-f(t j )At j )ij>(t j ), 
so that the Evolution Operator might be the strong operator limit 

U(t",t>)= Km 11(1-/(1^):= J] [1-/(0*]. 

1 l_ * j t">t>f 

If f(t) is a standard pseudo differential operator then the partial 
product approximations are pseudodifferential operators again. How- 
ever, if the order of fit) is positive, then the order of these pseudodiffer- 
ential operator approximations increases to infinity and the convergence 
of their symbols is out of the control. 

Fortunately, the backward Euler approximation 
suggests 

u(t",t')= n (l+fmr 1 

t">t>t> 

with zero order approximation symbols. 

Our main result is that for apt functions / this backward approxi- 
mation entails (in the spirit of the DFT ansatz) 



u(t"X)= n [(i - /(Odfrr 

f>t>t> 

leading to a Path Integral representation of the symbol <p\U(t", t')\q>. 

Incidentally, the Green function (the coordinate propagator) can be 
easily expressed via the symbol: 

<q"\U(t",t')\q'>= J dp< q "\p><p\U(t",t')\ q '>. 
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2 Rigorized Q-symbolic calculus 



This section provides necessary technical tools. 

For z = (q,p) G M. d x M. d introduce the complex coordinates z + = 
2~ 1 / 2 (q + ip), z~ = 2~ l l 2 (q — ip) so that the standard symplectic form 
[ipi, Qi), (P2, q2)} = V\(h ~ Viq\ on R 2d becomes 



Heisenberg Canonical Commutation Relations z \— > z over R 2d 
in a Hilbert space H is a linear map of M 2d to essentially self-adjoint op- 
erators on a common invariant subspace Q of TC (the Garding domain) 
such that [zi,z 2 ] = ftfzi, z 2 ]l. 

E.g., for the Schrodinger (position) representation on C 2 (M, d ) 



the Schwartz space (S(lR d ) may be chosen as a Garding domain. 

Other examples are the momentum or mixed momentum-position 
representations, holomorphic Bargmann-Segal representation (conducive 
to the coherent states Path Integral), Gelfand-Zak representation in the 
Solid State Physics, the Cartier compact representation in ^-functions. 

By the von Neumann-Stone theorem, for given h > any Heisenberg 
Canonical Commutation Relations is unitary equivalent to a direct sum 
of the Schrodinger representations. Thus we may chose the Garding 
domain Q^H) to be unitary equivalent to a direct sum of the spaces 



Correspondingly, the dual space Q'(TC) is unitary equivalent to a 
direct sum of the spaces <S'(IR d ). 



t[zi, z 2 ] — -{z\z 2 — z l z%). 



The h -Symplectic Fourier transform is defined as 




zip(x) = (qx)ip(x) + 



hdip 
i dx 




S{R d ). 
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Weyl operators / on 7i associated with generalized functions / G 
M. 2d are continuous linear operators from Q(7i) to Q'(TC) 

f = J /(C)expi[C,l]dA ft (C) 

wherein 

A version of the Schwartz Kernel Theorem states that a linear operator 
from Q{K) to Q'(TC) is continuous if and only if it is a Weyl operator /. 

Q -symbols. Consider a formal power series over C 

\a\>0 

A formal Q-symbol of / e S'(M. 2d ) is the formal power series over iS'(!R 2d ) 
defined via 

f n (C) = //0(C). 

Obviously, this has sense for polynomial f(z) when various Q provide 
common ordering rules according to the following table (c.f. [|TJ): 
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Name 


n(c) 


Ordering (d — 1) 


Weyl 


1 


q n p m <-» 


Standard 
(qp or Kohn-Nirenberg) 


e |[(c+) 2 -(r) 2 ) 




Antistandard 
(or pq) 


e -|[(C+) 2 -(C-) 2 )] 




Normal 
(or Wick) 




(^ + ) n (^-) m <-» 
(5 + ) n (i-) m 


Antinormal 
(or Anti-Wick) 




(2 + ) n (^) m <-» 
(z ) m (z + ) n 


Symmetric 


cosi[(C + ) 2 -(D 2 ] 


q n p m <-» 
\{q n p m +p m q n ) 


Born- Jordan 


sini[(C+) 2 -(C-) 2 ] 

|[(C+) 2 -(C-) 2 ] 


q n p m <— > 



7 



Suppose now that Q{() + for all C e K 2d . Then /(C)AXC) is 
meaningful for / G 5'(1R M ) if and only if 1/Q is a multiplier in iS(IR 2ci ). 
In such / n is called the strict fl-symbol of the distibution / G 

iS'(!R 2d ). E.g. every / G iS'(IR 2d ) has strict standard , antistandard and 
normal symbols. More generally f Q is called the strict Q -symbol of / 
if only n(C)f n (C) e 5'(M 2a! ). 

Of course not every Weyl operator has either antinormal, or sym- 
metric, or Born- Jordan symbol. 

Quasi-polynomials. Define (c.f. [HJ, Appendix 2) for m = (m 1; m 2 ), r = 
(ri,r 2 ), ri > 0, r 2 < 1/2, the class S(m,r) of quasi-polynomial f = 
{{fa{z), 0<h< h(f)} in such that 

d a J = Q (1)(1 + ^Dmi-nM/ima-raM 

wherein <9 Z = d/dz and a is a multiindex. 

As usual, S(— oo) = flS'(m,r), 5(oo) = U5'(m, r). 

A quasi-polynomial / is said to be asymptotic to a series 5^ a>At / Q 

a>fi 

with / Q G S(m a ,r a ), m a \ — oo, r Q \ — oo if for all v 

f -Y,F£S{m v ,r v ). 

The classical Borel-Hormander construction leads to the following 

Proposition. For every f G S(m, r) and a formal Q there is a g G 
S(m,r) asymptotic to [l/£l(^d z )]f. 

Such function g is called an asymptotic symbol f n of /. It is defined 
mod S(—oo). 

fi-products of quasi-polynomials. If fj are quasi-polynomials then 
fj act from Q(TC) to Q{7i) and therefore from <5'{7i) to Q'(TC), so that 
/ = /i/ 2 . . . /at is well defined. Actually / is quasi-polynomial, and 

TV 

J 3=1 
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wherein 

/C S2 (Ci, . . . , Cn) - c ) exp j - }JC„ a] > , 

(=, if f2 is strict). Generally, the integral is distributional, but abso- 
lutely converges for some Q, e.g., the normal one. 

The integral representation entails the asymptotic expansion 

wherein = + 77^), = 7^(7^ - 77^)- 

Trace. A density operator p : Q' — > is a Weyl operator with 
p G 5(IR 2a! ). The operator trace of / Q p n is well defined and may be 
evaluated for strict Q via the Trace formula: Tr(f n p n ) =< f n \p n > . 

3 Main Theorem. 

A quasi-polynomial / G S(m, r), m > 0, is called apt if for sufficiently 
small h it satisfies the following three conditions uniformly: 

• Quasi- dissipativity: Re(if) > 5, a constant . 

• Hypoellipticity: for all multi- indices a and < t', t" < T 

d a J{t",z) = O a (l)\if(t')-5\(l + \z\)- r ^h- r ^. 



• t-Continuity of f(t,- ) in S(m,r). 

Law of Inertia: If / is apt then all its asymptotic symbols f n are apt 
as well, albeit on different intervals of h. 

Also if / is hypoelliptic and real then f(t) are essentially self-adoint 
(c.f. [|19|, Proposition A2.1) in H. 
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Main Theorem. If f is an apt quasi-polynomial, then for suffi- 
ciently small ft, 

(1) The Cauchy Problem 

^ + /(t)V>0M) = O, ^(*,O) = <0o, 0<t<T, 

is proper on 7i relative G(H). 

(2) The evolution operator is the strong product integral 

u(t",f)= n [(l+^onp. 

t">t>t> 

(3) A strict Q -symbol u n (t",t',z) of the evolution operator is the 
limit in S'(M. 2d ) of the strict Q-symbols u v (t", t', z) of the partial 
operator products W t u >t >t i [(l + ^f(t J ,z))^)ras\V\^0. 

Proof (outline). We apply the f2-calculus along with the theory of 
Abstract Cauchy Problems (c.f. 0) and the theory of Finite Difference 
Methods for Initial Value Problems (c.f. [IB]) with the terminology 
thereof. 

The following statements hold for various intervals of positive ft. 
By the Law of Inertia, the anti-normal symbol of / is quasi-dissipative. 
Then (c.f. pf, Proposition 24.1) the real part of <il>\5 1 l + if(t)\ij» 
is greater than 7< ip\ip> with some constants 7 > and 8\. It is safe 



to assume that 8\ = 0. Together with the hypoellipticity (c.f. |^9f| , 
Theorem 25.4) this entails that ||[A1 + if(t)]~ 1 \\ < 1/A for positive A so 
that the operators f(t) (l,0)-stable family in 7i. 

When both ip and f{t)ip belong to 7i for some t = t then it is 
so for all t by the virtue of the hypoellipticity. The space T of all 
such i[) is dense in 7i and is a Hilbert space relative the new Hermitean 
product < V# >o=< f(to)i/j\f(t )ip >. Now the f(t) : T -* H form 
a t-continuous family of bounded operators. Moreover, < f(t)ip\ip>o 
=<g(t)ip\ip> with g(t)=f(t ) : ' i f(t )f(t) so that g is apt again and thus 
(as above) is (l,0)-stable in T . By the Hille-Yosida theorem flQ| , f(t) 
generates for every t a contractive operator semi-group in T . Since 
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the family f(t) — f(T) has similar properties, the theorem 7.7.13 of || 
establishes (1), the properness of the Cauchy problem. 

This leads to a preliminary Product Integral representation U (t", t') = 

n t »> t >t' C 1 + x/" 1 ] ( ci - the P roof of the theorem 7.7.5 of §). It im- 
plies the Product Integral representation (2) via the Lax Equivalence 
Theorem [16] whereby the required consistency is checked via the Weyl 
calculus. 

The last statement (3) follows from the trace formula. 



4 Path Integrals in Quantum Field The- 
ory. 

Infinite dimensional phase spaces. In the case of d = oo there 
are non-isomorphic phase spaces and the symplectic structures usually 
appear with extra features. 

Our phase space is based on a separable Frechet nuclear space Z 
over C with a "dotless" hermitian product zw*. 

If TC is the corresponding Hilbert space completion of the Z, and 
Z* is the corresponding anti-dual of Z then 

Z ^ H ^ Z* 

is a Gelfand nuclear triplet. 

The phase space is Z taken over K with the symplectic form —lm(zw*). 
It is also a pre-Hilbert space with the scalar product Ke(zw*). 

Complex Gaussian rigging. (C.f. ||.) The Gaussian measure 7^ 
of covariance 1/h is defined via its characteristic function 

f e lRe{zC) d lh (C) = e~ zz * /2h } zeZ, 
Jz* 

so it stands for the non-existent (^r)°° exp(-h(*(*)d(* . 

The Bargmann-Segal space [H] is the closure of the subspace of the 
continuous complex analytic polynomials on Z* in C 2 (Z*). Its elements 
are entire functions h(z*) of order 2 and type < h/2: 

h(z*) = 0(e hp ^ )2/2 ) 
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for some dual semi- norm p on Z*. 

Let [Z] denote the space of entire functions h(z*) of order 2 and 
minimal type, and [Z\* denote the space of entire functions h(z*) of 
order 2 and maximal type. Then [Z] is naturally a separable nuclear 
Frechet space, and [Z\* its anti-dual. Thus 

[Z] - [H] - [2]* 

is another Gelfand triplet, the Complex Gaussian rigging of the triplet 

z ^ n ^ z*. 

The coherent states e w (z*) := exp(—hwz*), w G Z, form a total (over- 
complete) set in [Z]. 

Free bosonic field over Z in [Z]*. Let z — > z be an antilinear 
conjugation on Z and correspondingly on Z*. Set 

^+ = z/V2 e z, z~ = z/V2 e z. 

The operators z + and z~ are defined on f £ [Z\* as 

5 + /(C) = (zC)f(C), r = hd z f(C). 

They represent the Canonical Commutation Relations (CCR): 

[z-,z + ] = hi. 

The coherent states are entire vectors for the CCR. 

Wick Operators are, by definition, the continuous linear operators 
W from [G] to [£]*. The Wick symbol of W is 

w(z + ,z~):=e- z+z - [ [We z+ (C)]e z -(C + )dln(0, z G Z. 
J z* 

The Wick symbol w is an entire function on Z x Z, so that the 
operator w := w(— z + ,z~) is well defined on the coherent states and 
W = w and W is its closure. 
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Q -symbols. Consider a formal complex power series 1 + J2\ a \>o za 
on Z*. The formal fi-symbol of w(() is 



w 



i/n(?0, 

2 



W (C). 



Quasi-polynomial w(() is the family {wa(C) ■ < h < h(w)} such 
that for a dual semi-norm p on Z* 

<9>(C) = +p(C)) mi_ri|o| ^ m2 " r2H . 

The class of such families is denoted S(m,r), m = (m 1 ,m 2 ),r = 
(ri,r 2 ). 



Weyl symbols and operators. The Weyl symbols f(z) of w(z) 

i, 

' 2- 



correspond to Q(z) = exp(— \z + z ) (c.f.the table above): 



f(z) = exp 

Conversely 



h 2 d 2 



W(z)= / f(z-()d lh y 2 ((). 

J z* 

(Note: not every Wick operator has a strict Weyl symbol.) 

The corresponding Wick operators are the Weyl operators f. The 
Borel-Hormander constructions for a countable fundamental family of 
dual gaussian semi-norms followed by the Cantor diagonal trick imply 
that for every Q and Weyl / G S(m, r) there is g G S(m, r) asymptotic 
to 

The fi-symbols of the operator product f\f 2 ---fN have the asymp- 
totic expansions just as in the case d < oo. However their integral 
representations are known rarely. Fortunately, for the normal symbols 

W of W\W 2 ...Wn 

/N N-l 
Y[ e z i z t-i Wj (zJ ,z+_ t ) Yl dihizjizf^), z£ := z + , z^ := z 

3=1 J=l 
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Finally, as in the finite-dimensional case, the Main Theorem holds 
in infinite dimensions (with the same proof) at least for the strict Wick 
symbols of evolution operators in G = \H}. In the latter case the sym- 
bol approximations are absolutely convergent multiple integrals with 
respect to d^n over the infinite-dimensional phase space Z. 



5 Conclusion and outlook. 

1. The phase space Path Integral (according to L.Shulman, "a diffi- 
cult form" of the path integral) was originated in different ways 
by Feynman himself J7[ in 1951 and by Tobocman in 1956. 
The Coherent State discretization was introduced in 1960 by 
Klauder in the Schrodinger representation and in 1962 by 



Schweber |L7| in tne Bargmann-Fock represenatation. In the 
70 's Berezin @] considered various discretizations on the basis 
of pseudo-differential analysis. However no convergence of the 
discretizations has been proved until now. 

On the other hand Daubeshies & Klauder []5| have established 
in 1984 that a wide class of coherent state path integrals (essen- 
tially with self-adjoint polynomial hamiltonians) on a flat finite- 
dimensional phase space may are limits of Wiener Integrals on the 
space of paths in the phase space. They even suggested that Feyn- 
man type time-slicing construction is impossible for the phase 
space Path Integrals. 

We have presented a rigorous time-slicing Phase Space Path Inte- 
gral construction for the symbols of the Evolution Operators with 
wide variety of smooth hamiltonians both in finite and infinite de- 
grees of freedom. The convergence is established only for small 
H, in agreement with the postulated semi-classical nature of the 
Path Integral which relates the classical and quantum dynamics. 

According to the f2-calculus, the discretizations of the Path In- 
tegral are distributional multiple integrals. E.g., as mentioned in 
the Introduction, the traditional discretization of the Phase Space 
Integral comes from the standard f2-calculus. Similarly, the Co- 
herent State Path Integral discretization is associated with the 
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normal fi-symbol in which case the multiple integrals are abso- 
lutely convergent. 



4. The last statement in the Main Theorem is equivalent to a mod- 
ified DTF-ansatz: the f2-symbol of the short time propagator is 
approximately equal to [1 + j^f(t, z)At}~ 1 . However, in the case 
of the normal f2-calculus (because of the absolute convergence) 
one may consistently replace it with the more customary ansatz 
exp[-if(t,q,p)At}. 

5. Our Path Integrals are "pathless", in agreement with the Uncer- 
tainty Principle: no quantum path in the phase space (c.f. |T3 



for an illuminating somewhat different point of view). Yet they 
are semiclassical in the following sense: the principal terms in the 
^-expansions of the partial products symbols are the backward 
Euler approximations of the corresponding classical Hamilton- 
Jacobi equations. 

6. We have rigorized the f2-calculus of Agarval & Wolf |l[ to jus- 
tify numerous Path Integral discretizations and as an important 
techniques. However in the infinite degrees we have been able to 
prove the convergence only for normal (Wick) symbols. 

Actually the formal fi-calculus is a special case of the formal *- 
calculus U . Since on the finite-dimensional flat symplectic space 
all ^-products are formally equivalent, our results yield a con- 
struction of the formal *- exponential, a solution of a well known 
problem (c.f [|18| for an interpretation of the Evolution Operator 



symbol clS cL ^-exponential). 

7. Most of the other mathematical interpretations of the Path In- 
tegral are primarily in terms of various distributional integrals 
on the paths in the configuration space : first, by Kac [JH 



via 



analytic continuation to a Wiener integral (the Feynman-Kac for- 
mula), followed by DeWitt-Morette [frj| in terms of prodistribu- 



tions, by Albeverio and H0eg-Krohn || in terms of the Parseval 
equation for the oscillatory Gaussian integrals, and by Hida & 
Streit || in terms of White Noise disributions. Notably, these 
Path Integrals are associated only with the Schrodinger hamilto- 
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nians (essentially) of quadratic growth, with the presumed "Feyn- 
man measure" built from the kinetic energy term. 
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